Spatial diffusion in a periodic optical lattice: revisiting the Sisyphus effect 
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We numerically study the spatial diffusion of an atomic 
cloud experiencing Sisyphus cooling in a three-dimensional 
lin_Llin optical lattice in a broad range of lattice parameters. 
In particular, we investigate the dependence on the size of the 
lattice sites which changes with the angle between the laser 
beams. We show that the steady-state temperature is largely 
independent of the lattice angle, but that the spatial diffusion 
changes significantly It is shown that the numerical results 
fulfil the Einstein relations of Brownian motion in the jumping 
regime as well as in the oscillating regime. We finally derive an 
effective Brownian motion model from first principles which 
gives good agreement with the simulations. 



model of Brownian motion [|2||14|]. In this regime, spa- 
tially averaged friction coefficients cto and momentum 
diffusion coefficients D p have been derived and the va- 
lidity of the Einstein relations 



k B T = — 

a 



and 



D„ = 



k B T 
ceo 



(1) 



(2) 



PACS. 32.80.Pj Optical cooling of atoms, trapping - 42.50.Vk 
Mechanical effects of light on atoms, molecules, electrons and 



I. INTRODUCTION 

Laser cooling and trapping was one of the major ad- 
vances of the last part of the 20th century. In 1997, 
the Nobel prize in physics was awarded to Steven Chu, 
Claude Cohen- Tannoudji and William D. Phillips for 
their works in this domain [jjj, in particular, for their 
discovery of the Sisyphus cooling effect j^j which permits 
to achieve sub-Doppler temperatures and is widely used 
in various different laser cooling schemes The Sisy- 
phus effect in optical lattices was studied in a large 
variety of systems and of field configurations. A large 
number of results were obtained on temperature, local- 
ization and spatial order Q and an excellent agreement 
between the experimental observations and the theoret- 
ical predictions was found. Much less work has been 
done to study the spatial diffusion in optical lattices, but 
also for this problem a reasonable agreement was found 
between the models and the experiments Jfl|,p|-pd| . 

However, no detailed study of the dependence of spa- 
tial diffusion on, for example, the different directions of 
an anisotropic lattice or on the size and shape of lattice 
sites has been performed so far. Very recently, optical 
lattices and atomic transport therein has attracted new 
attention with the study of quantum chaos |T^ | and the 
achievement of Bose-Einstein condensation by purely op- 



where T is the steady-state temperature, k B the Boltz- 
mann constant, and D s the spatial diffusion coefficient, 
has been shown. On the contrary, in the so-called oscil- 
lating regime, where an atom travels over several optical 
wavelengths before being optically pumped into another 
internal state, the friction force has been shown to be 
velocity-dependent , 



F(v) = 



-a v 



i + (v/v c y 



(3) 



tical means |13 



Spatial diffusion in one-dimensional (ID) Sisyphus 
cooling schemes is fairly well understood. In the so-called 
jumping regime, where an atom undergoes several opti- 
cal pumping cycles while moving over one optical wave- 
length, the atomic motion can be understood by a simple 



where v c denotes the capture velocity of Sisyphus cool- 
ing. In this situation, an analytical derivation of the 
spatial diffusion coefficient has still be found ||, but an 
interpretation in terms of a simple Brownian motion no 
longer works. In particular, it is found in Ref. || that 
the behavior of the spatial diffusion coefficient as a func- 
tion of the atom-light interaction parameters (laser light 
intensity and detuning) is dramatically different in the 
oscillating regime compared to the jumping regime. 

In higher dimensional setups, such a difference of the 
spatial diffusion behaviors in the jumping and the oscil- 
lating regimes is expected, too. Moreover, the important 
difference of the mean free path of a diffusing atom in 
these regimes may induce different behaviors of the spa- 
tial diffusion coefficients as a function of the lattice peri- 
ods. These are the main issues considered in the present 
paper. 

In this work we perform a detailed study of spatial dif- 
fusion in the so-called 3D-lin_Llin lattice [jl5| . Using semi- 
classical Monte-Carlo simulations wc find that equalities 
of the form of Eqs. ([!]) and (|^) still hold in the oscillating 
regime. This suggests that an interpretation by a Brow- 
nian motion should still be possible. We derive such a 
model from basic principles assuming a thermal spatial 
distribution and taking into account some specific prop- 
erties of our optical lattice and find a good quantitative 
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agreement with the numerical results. In particular, we 
calculate an effective friction coefficient a in a range of 
parameters containing both the jumping and the oscil- 
lating regimes. We find that in the oscillating regime, a 
increases with the lattice beam intensity and decreases 
when |A| increases, in strong opposition to the friction 
coefficient ao calculated in the jumping regime. 

Our work is organized as follows. In Sec. || we de- 
scribe the specific laser and atom configuration for the 
3D optical lattice that we consider here and discuss sev- 
eral important features of the optical potential surfaces. 
In Sec. [II, we present a physical picture of spatial diffu- 
sion in periodic optical lattices and we particularly fore- 
cast a dramatic change of the behavior of the spatial 
diffusion coefficients not only versus the atom-light in- 
teraction parameters but also versus the geometrical pa- 
rameters (spatial lattice periods) when going from the 
jumping to the oscillating regime. In Sec. IV we derive 



an effective Brownian motion model which we compare 
in the following sections with the numerical results on 
the steady-state temperature (Sec. |v|) and on the spatial 
diffusion (Sec. VI). Numerical res ults on the friction co- 
efficient are then discussed in Sec. |VI]| and the validity of 
the Einstein relations is shown. Finally, we summarize 



our results in Sec. VIII 



II. 



SISYPHUS COOLING IN THE 3D-LIN±LIN 
CONFIGURATION 



so 



a 2 + r 2 /4 



«1, 



(4) 



where so is the saturation parameter for an atomic tran- 
sition with a Clebsch-Gordan coefficient of one, f^o is 
the Rabi frequency for one laser beam, A the detuning 
of the laser beams from the atomic resonance frequency, 
and r the natural width of the atomic excited state. This 
domain is known to lead to the lowest temperatures. In 
this situation we may adiabatically eliminate the Zeeman 
sublevels of the excited state, leading to a theory which 
only involves the ground state sublevels |±) of angular 
momentum ±1/2 [Q . 




The Sisyphus effect cools a cloud of multi-level atoms 
when a laser field induces spatially modulated optical po- 
tentials and pumping rates in such a way that a moving 
atom on average climbs up potential hills before it is op- 
tically pumped into a lower lying potential surface ^ . In 
this case kinetic energy is converted into potential energy 
which is subsequently carried away by a spontaneously 
emitted photon, thereby reducing the total atomic en- 
ergy. 

In this paper we study the so-called 3D-lin_Llin config- 
uration |lq| . It is obtained from the standard lD-lin_Llin 
configuration Q by symmetrically splitting each of the 
two laser beams into two parts at an angle 9 X and 9 y , 
respectively, with the (Oz) axis in the (Oxz) and (Oyz) 
planes respectively. The resulting configuration consists 
of two pairs of laser beams in the (Oxz) plane and in the 
(Oyz) plane, respectively, as depicted in Fig. |], with or- 
thogonal linear polarizations. An important property of 
this configuration in contrast to 3D setups built of more 
than four laser beams is that the interference pattern 
and thus the topography of the lattice does not change 
because of fluctuations of the relative phase between the 
various laser beams. Instead, such fluctuations only in- 
duce displacements of the lattice. 

As in most theoretical work, we will consider atoms 
with a ground state of angular momentum 1/2 and an 
excited state of angular momentum 3/2. Experiments 
usually fall into the low saturation regime defined by 



FIG. 1. Laser beam configuration for the three-dimensional 
lin_Llin Sisyphus cooling. Two pairs of laser beams with 
crossed linear polarizations induce polarization and light shift 
gradients in the three directions of space. 

An atom in state |±) then experiences an optical po- 
tential U± given by 



U±(x,y,z) 



8hA' 



cos(k x x) 2 + cos(k y yY 



=F cos^a;) cos(k y y) cos(k z z) 



where 



A = A^ 
2 



is the light shift per beam and 



k x = ksin(9 x ), 

ky = k Sill(9y), 

k z = k[cos (9 X ) + cos(6> y )], 



(5) 



(6) 



(7a) 
(7b) 
(7c) 



with k the laser wavenumber. The optical potentials are 
then periodic in the three directions of space with periods 
A,; = 2n/ki. Equation (|^) shows that U±(x 1 y,z) has the 
same functional dependence on k x x and on k y y but a 
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different one on k z z. We will therefore in the remainder 
of the paper concentrate on a two-dimensional subsystem 
only depending on x and z while fixing y = 0. Previous 
comparisons between ID and 2D models have shown that 
the general behavior of the dynamic variables is the same 
in different dimensions but that scaling factors appear 
||. We thus expect that the results of our work give 
physical interpretations to full 3D laser cooling schemes 
but exact numerical values will be changed. Moreover, we 
will assume a single lattice angle 9 — 9 X — 8 y since this 
gives rise to a vanishing mean radiation pressure force in 
all directions. The general shape of the optical potential 
is plotted in Fig. 0. 



because of the geometrical dependence of the border be- 
tween the jumping and the oscillating regimes, 



r 



A, 
A 



(10) 



The asymmetry between the x and z directions can 
be seen most easily in a plot of the optical potentials 
along the x and z axis, respectively, as shown in Fig. ||. 
These have different shapes and particularly the crossing 
between both potential curves is higher in the transverse 
direction (x) than it is in the longitudinal one (z). As we 
will see later, this induces significant differences in the 
cooling and diffusion properties. 



(X,Z) ffl\ A'o I 




FIG. 2. Section at y=0 of the optical potential associated 
with the Zeeman sublevel |— ) of the atomic internal ground 
state in the case of a negative detuning (A < 0). U+ can 
be obtained from U- by a translation of half a spatial period 
along one of the axes (Ox) and (Oz). 

In the 3D-lin_Llin configuration, the bottom of each 
potential well is harmonic in first approximation with 
main axis x, y and z and with the following frequencies: 



n XjW = 4^/|A&|w r -|* 



(8a) 
(8b) 



where lor = hk 2 /(2M) denotes the recoil frequency. The 
optical pumping time is 



9 

8TV 



(9) 



where T' Q = 



Tsq/2 is the optical pumping rate. The 
jumping regime corresponds to a domain where Qi<til/T p , 
that is, when an atom undergoes many optical pumping 
cycles during a single oscillation or during a flight over 
a single potential well. On the contrary, the oscillating 
regime corresponds to r2i3>l/r p . In this case, an atom 
can oscillate or travel over many wells without undergo- 
ing any pumping cycle. Note that in a 3D-lin_Llin optical 
lattice, the regimes can be different in different directions 



U± (x,y=0,z=0)/fi|A'o 




k x x 



k z z 



FIG. 3. Sections at (a) y 
optical potential surfaces U-\ 
curve) . 



= z — 0, (b) x = y = of the 
(solid curve) and U- (dashed 



As the starting point for the theory presented here we 
use the standard Fokker-Planck equation (FPE) of the 
semiclassical laser cooling theory ^|l7],[l8| where the ex- 
ternal degrees of freedom of the atoms are treated as 
classical variables. This is obtained from the Wigner 



transform |19| of the full quantum master equation for 
external as well as internal degrees of freedom under the 
assumption of a momentum distribution which is much 
broader than a single photon momentum, AP 3> fik. 
The FPE for the populations II±(r,p) of |±) reads 



Pi 



dt + ^di - (diU±)d Pi 



M 

7 T n T - 7 ± n ± 

-Fi ± d Pi U ± 



n± 



+D l i ± d Pi d Pj n ± + D^ ± a Pi d Pj u T . 



(ii) 



Here i, j — x,z and summation over i and j is assumed. 
In this equation, j± is the jumping rate from the Zeeman 



3 



sublevel |±) to the sublevel |=f), F± ± represents the ra- 
diation pressure force and D % ± ± the momentum diffusion 

matrix for atoms in the internal state |±). F±zc. an d 
are the corresponding quantities associated with jumps 
between different internal states 16 1. Note that all these 
coefficients only depend on the atomic spatial position 

A numerical solution of the FPE can be obtained by 
averaging over many realizations of the corresponding 
Lan-gevin equations pc[ . Such a realization consists in 
following the trajectory of a single atom which jumps 
between the two optical potential surfaces corresponding 
to the two internal states with the appropriate probabil- 
ities. Between subsequent jumps the atom experiences 
potential and radiation pressure forces as well as random 
momentum kicks which mimic the momentum diffusion 
according to the coefficients of the FPE. We have per- 
formed a large number of such semiclassical Monte-Carlo 
simulations in order to investigate the dependence of the 
steady-state temperature, the friction coefficient, and in 
particular the spatial diffusion coefficient on the various 
lattice parameters such as detuning A, light shift A , and 
lattice angle 6. We will disc uss these numerical results 
later in Sees. [v|, [v| andfvTi. 



III. 



PHYSICAL PICTURE OF SPATIAL 
DIFFUSION 



As a result of the Sisyphus effect, the atoms are cooled 
and trapped in the potential wells and optical lattices are 
usually described as atoms well confined in regularly ar- 
ranged sites (note that in the lin_Llin lattice, the spatial 
periods of these sites are A,/2, the trapping sites cor- 
responding alternatively to U+ and U- potential wells). 
However, in bright optical lattices the atom confinement 
is not perfect because of the strong interaction with the 
laser light. Two different processes then produce atomic 
displacements between different trapping sites, inducing 
spatial diffusion (see Fig. ^). For the sake of simplicity, 
we describe these processes in one dimension but they 
occur analogously in higher dimensional setups. 

On the one hand (see Fig. || left), a trapped atom still 
undergoes fluorescence cycles and thus takes random re- 
coils due to photon absorption and re-emission. Hence, 
the oscillating motion of the atom gets perturbed. Par- 
ticularly, the atom can explore regions where its poten- 
tial energy (U± if the atomic internal state is |±)) is not 
minimum (U± > U^). In such regions, optical pumping 
cycles preferentially transfer the atom into the lower po- 
tential curve and the atom is cooled and trapped in the 
neighboring potential well (elementary Sisyphus cooling 
process). This process induces atomic transfers from a 
site to a neighboring one in another potential curve. 

On the other hand (see Fig [3| right), when a trapped 
atom oscillates in a potential well, it has a small but non 
zero probability of being optically pumped into the up- 



per potential curve. In the jumping regime the atom is 
immediately pumped back into its initial trapping po- 
tential well. This effect thus induces heating and noisy 
oscillating trajectory of the atom leading indirectly to 
spatial diffusion via transfers between neighbouring po- 
tential wells. On the contrary, in the oscillating regime 
the atom is not immediately re-pumped and travels over 
several potential wells before undergoing an elementary 
Sisyphus cooling process again which traps it into an- 
other potential well. 




FIG. 4. Processes of spatial diffusion in an optical lattice. 
Left: process due to the heating. Right: process due to optical 
pumping. 

The diffusion process linked to optical pumping is 
much more efficient than the one due to recoils except 
for very small laser detunings jw| . We will thus focus on 
the second process to describe spatial diffusion in peri- 
odic multi-dimensional optical lattices. We will see that 
the differences of this process in the jumping and the 
oscillating regimes induce a dramatic difference in the 
behavior of the spatial diffusion coefficients. 

In a simple model, we can describe the diffusive behav- 
ior of the atomic cloud as random walks of atoms between 
periodic trapping sites pi]]. Let us assume that an atom 
is trapped in one specific potential well and jumps after 
a time r to another well. The spatial diffusion coefficient 
in the direction i £ {x, y, z} is 



Dsi = — 
b 2r 



(12) 



where di is the mean free path in direction i. 

In the jumping regime, as discussed above, an atom 
essentially transfers from a trapping site to a neighboring 
one and thus d{ ~ A^. The life time of an oscillatory 
external state is on average of the order of l/T' [ p2[ , 
independently of A,;. Hence, r ~ 1/Tq and 



D 



Si 



A?r'. 



(13) 



In the oscillating regime, an atom travels over several 
lattice sites before it is trapped again. Here d ~ vt 
where v ~ \Jkj^TjM is the average velocity with M the 
atomic mass and r the time of one flight. As it will be 
justified in Sees. IV and |v|, UbT is proportional to 7j|A | 
independently of A and 9. r is of the order of 1/T' again 
p2|. Hence, 



(14) 



4 



Note that we implicitly assume straight line flights and 
we do not consider anisotropic effects. Obviously, be- 
cause of the potential shape and anisotropy, a dependence 
of Dsi on X x , X y and X z should be added in Eq. (|14) 
but our simple model does not provide its determina- 
tion. Nevertheless, in the case of a strong anisotropy 
(A, <C Xj), the setup is almost one-dimensional in the 
i-direction so it is expected that D$i does not depend 
on Xi. Indeed, at large space scale, the length of one 
flight of a particle moving on a ID-periodic potential is 
independent of the periodicity. 

We want to emphasize that these discussions are only 
valid for lattice parameters far away from the domain 
of decrochage for spatial diffusion. For small potential 
depths, rare long flights dominate the diffusion which 
therefore becomes anomalous 0. 



IV. BROWNIAN MOTION MODEL 

Before turning to a detailed discussion of the numeri- 
cal results, we will now derive a simple analytical model 
which will help to understand the main features of the 
cooling scheme. 

The basic idea of this model is to consider the atomic 
dynamics as a Brownian motion, as it has been success- 
fully applied to Doppler cooling [H and ID Sisyphus 
cooling in the jumping regime [0Jlj], but taking some 
properties of localization into account |^3|,^4| . The fun- 
damental ingredients are therefore the derivation of an 
average friction force and an average momentum diffu- 
sion coefficient. 

Let us consider an atom at position r = xe x + ze z 
moving with a constant velocity v = v x e x + v z e z . In this 
case the FPE (nil) reduces to 



{d t + vv)n± = - 7± n± + 7T n T 

= -(7- +7+)n± +7t 



(15) 



where we used IL 



n 



1. The jump rates are 



7± = gT^l + cos 2 k x x ± 2 cos k x x cos k z z). (16) 



Expanding the populations in powers of the velocity in 
the form II ± = ^ n II^ and inserting into Eq. (15), 
yields the stationary solutions 



IT 



(0) 



It 



n 



with the operator 



(n) 



7- +7+ 

Q™r4 0) 



7- +7+ 



rV. 



(17a) 
(17b) 

(18) 



Formally, the velocity and position dependent level pop- 
ulations can thus be written as 



n± 



l-O 



(0) 



(19) 



The total force averaged over the internal atomic states 
is then given by 



F(r,v) = -(IT+VC/+ +n_Vf/_). 



(20) 



In order to derive a space and velocity averaged friction 
force in the form F = — (v x a x e x + v z a z e z ) we now have 
to make certain assumptions on the stationary atomic 
distribution. 

In ID laser cooling one usually assumes a flat spatial 
distribution of the atoms in the lattice. However, as we 
have seen in Sec. [0], the shape of the 3D optical po- 
tential differs significantly in the different directions. In 
particular the potential barrier between neighboring po- 
tential wells is much higher in the transverse than in the 
longitudinal direction. From Eq. (^|) and Fig. ^ we see 
that in the x direction the potential depth is of the order 
of 167i| Aq|/3 whereas we will see later that the steady- 
state temperature is of the order of 2?i|Ag|. Therefore 
we expect strong localization in that direction. On the 
contrary, in the z direction the potential depth is of the 
order of 8?i|Aq|/3 and the atoms will be less localized. 
Instead of a flat spatial distribution of the atoms we will 
thus assume a thermal distribution 



P(r) oc exp 



1 -U+ + 1+ U- 
k B T( 1+ +7_) 



(21) 



corresponding to the optical potential averaged over the 
internal atomic state and for a given, yet unknown, tem- 
perature T. Let us emphasize that assuming a thermal 
distribution is not a priori justified for laser cooled sam- 
ples but show significant quantitative deviations from 
such a simple behavior |25| , |26|1 . In fact, our numerical 
simulations give actual spatial distributions in qualitative 
agreement with Eq. |l| but show significant quantitative 
deviations from such a simple behavior. However, we 
only use this approximation here to obtain a qualitative 
understanding of the exact results obtained numerically 
and we will see later that our results derived here are in 
good quantitative agreement with the simulations. 

Because of the symmetry of -P(r), only terms contain- 
ing odd powers of the velocity in H± contribute to the 
averaged force. We may thus restrict ourselves to 



j^odd 



En 



(2r, 



1 



1 



Q* n± 



1) 



(22) 



As a further simplification we will now also average over 
velocity and therefore replace Q 2 by 



(Q 2 



2ujB,kBT/h 
<(7-+7+) 2 > 



{{k x /kf + [k z /kf 



(23) 



where (...) is the spatial average with respect to -P(r). 
Equation (|23| ) is obtained by substituting vi — > 0, vf — > 
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k B T/M, df -> -kf (i = x,z) in Eq. (|T|). From this wc 
finally get the friction coefficients 



a, 



6 ~ fiki Ci 



with 



C x — 

c z = 



i - (Q 2 ) 

sin 4 (k x x) cos 2 (k z z) 
[l + cos 2 ^)] 3 /' 

cos 2 (fc x x) sin 2 (fc z z) ' 
[1 + cos 2 {k x x)} 2 



(24) 

(25a) 
(25b) 



Along the lines of [f27f and using the same approxima- 
tions as for the friction, we derive averaged momentum 
diffusion coefficients 



_ 8±A' h 2 k 2 C< 
Pl ~ 1 - (Q 2 ) 



with 



C = ( in (0) n (0) s[n2 ( kxX "> cos2 ( fc ^) 

+ 1 + cos 2 (k x x) 

C = ( 4n (0) n (0) cos2 ( kxX ^> sin2 ( fc ^) 

+ 1 + cos 2 (k x x) 



(26) 

(27a) 
(27b) 



For simplicity we have not written a term of the mo- 
mentum diffusion which arises from the random recoil 
of absorbed and spontaneously emitted photons. This 
varies as T' and thus can be neglected for large detuning 

(|A|>r). 

In a model of Brownian motion, the steady-state tem- 
perature then fulfills 



k B T = - (Be. + Bi 

2 \ «i a z 



(28) 



that is, averaging Eq. (|l|) over the x and z directions. 
The right hand side is obtained from Eqs. (|24|) and 
( p6| ) which themselves depend on the temperature via 
Eq. (|2l]). Thus, Eq. ( |2§| ) yields an implicit equation for 
kBT/h\A' Q \ which can be solved numerically, e.g., by it- 
eration of Eqs. (pl"l)-(^8|) recursively. Note that no lattice 
parameter, such as the lattice angle or the laser detuning, 
appears in this equation. T is thus strictly proportional 
to ?i|Ag| and independent of 9 and A. We find 



k B T = 1.545 fi\A' \ 



(29) 



This temperature can then be used to determine the fol- 
lowing numerical values: 



0.0356, 
0.0409, 



C z 



0.0747, 
0.0874, 



and 



(( 7 -+7+) 2 ) =2.55r^. 



(30) 



(31) 



Finally, we want to derive an approximate expression 
for the spatial diffusion coefficients. To this end, we 
must again take the atomic localization into account. 
While Eq. (0) for the relation between temperature, mo- 
mentum diffusion coefficient and friction coefficient ap- 
proximately holds for trapped and untrapped atoms, 
the corresponding equation (j^) for the spatial diffusion 
only holds for free atoms. Indeed a cloud of completely 
trapped atoms achieves a stationary spatial distribution 
and hence shows no spatial diffusion. Using Eq. ( |29| ) 
and the assumption of a thermal momentum and spatial 
distribution, we calculate numerically that a fraction of 
55.6% of all atoms have a total energy above the potential 
depth along z, and a fraction of 15.3% above the poten- 
tial depth along x as discussed before. Taking only these 
free atoms into account, we finally obtain effective spatial 
diffusion coefficients, which correspond to those observed 
in the numerical simulations or in actual experiments, in 
the form 



D Sx = 0.554— I — 



M \k x 

uj r a r 



D 



Sz 



0.956 



M \ k 



^JK\ + 1 - 2l ^ [{kz/k)2 + {K/k)2] > ' ,:!2) 



This expression is in good qualitative agreement with our 
physical discussion (see Sec. HI) in both the jumping and 
the oscillating regime. We will further discuss this point 
in Sec. [VI. 



V. STEADY-STATE KINETIC TEMPERATURE 

We performed a systematic study of the temperature 
and the spatial diffusion as a function of the lattice pa- 
rameters, exploring a large domain containing both the 
jumping and the oscillating regime. More precisely, we 
performed numerical semi-classical Monte-Carlo simula- 
tions with the following parameters: 

6 = 15°, 30°, 45°, 60°, 75°, 

A/r = -2, -3, -5, -10, -15, -20, -25, -30, 

Ao/u; r = -150, -300, -450, -600, -750. 

In this and the following sections we present the results 
of the simulations and compare them with the analyti- 
cal model discussed above. All of our discussions and 
conclusions rely on the complete set of data, even if the 
figures only contain a few sample curves for the sake of 
clarity. 

We first study the atomic cloud steady-state tempera- 
ture resulting from the competition between slowing and 
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heating processes. We calculate the average square ve- 
locity over the whole atomic cloud at each time step of 
the simulations. The temperature in direction i is then 
given by 



k B T % = M(i 



(33) 



A second method to obtain the temperature numerically 
is to fit a Gaussian to the simulated momentum distribu- 
tion. We have checked that the widths of these Gaussians 
indeed give the same temperatures as those obtained 
from Eq. (33). For a broad initial velocity distribution, 
the temperature first decreases in time (thermalization 
phase), but finally reaches a steady state. 

The thermalization time lies between 50/Fq and 
100/Fq for large enough lattice angles (9 > 30° in the 
simulations), but for 9 = 15° it is approximately 400/rd 
in the transverse direction. In this case, the spatial pe- 
riod is large and an atom needs to fly a long time to 
undergo efficient Sisyphus cooling. Such an increase of 
the cooling time is also observed in the longitudinal di- 
rection for large angles but is less important because A z 
does not reach very large values. 




15 20 25 
A| /r 



FIG. 5. Transverse (a) and longitudinal (b) steady-state 
temperatures versus the laser detuning A for a fixed lat- 
tice angle 9 — 45° and for various light shifts (the values of 
|Ao|/uv ar e quoted to the right of the corresponding curves). 

Both the transverse and the longitudinal temperatures 
are shown in Fig. versus the laser detuning for a fixed 
lattice angle and for several light shifts per beam. This 
exhibits two domains where the temperature behaves dif- 
ferently. For small detunings (| A| < 10r) we find a rapid 
decrease of the temperature with increasing |A|, whereas 



for large detunings the temperature is independent of 
A but increases approximately linearly with Aq. This 
agrees well with the general form 



A, + Bi 



(34) 



where Ai, Bi and Cj (i — x,z) are numerical factors, 
as, for example, has been found in Rcf. j^]. The first 
term of Eq. (gj) has not been found in Sec. IV because 
in the momentum diffusion, Eq. (|26|), we neglected the 
term due to absorption and spontaneous emission. 

Let us first concentrate on the oscillating regime, 
|A| 3> r. In Fig. H we plot the temperature versus the 
lattice spatial period for various laser detunings and light 
shifts in this domain. We find that the temperature is 
nearly independent of the lattice angle, i.e., of the spa- 
tial periods. The temperature is thus strictly linear in 
the potential depth and independent of any other lattice 
parameter as predicted by Eq. (p9[ ) . Such a property has 
been observed experimentally Pql but we emphasize that 
our result here holds for a broader range of parameters in 
the oscillating regime as well as in the jumping regime. 
A linear fit to the numerical results gives 



k B T x ~ 2^|A^| + 123.37iu; r , 
k B T z ~ 1.437i|A£,| + 62.6ftw r . 



(35a) 
(35b) 
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FIG. 6. Transverse (a) and longitudinal (b) temperatures 
versus the corresponding lattice period for various laser de- 
tunings (triangles: A = — 10F, squares: A = — 20F, circles: 
A = — 30r) and light shifts (the values of lAdl/av are quoted 
to the right of the corresponding curves) . 
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In the range of parameters investigated here, these val- 
ues of the temperatures agree with Eq. (|2^) with an ac- 
curacy of about 10%. 

Equations (35) show that the temperature is 
anisotropic, T x > T z , which is in good agreement with 
experimental observations by A. Kastberg and coworkers 
p9| . This is a consequence of the asymmetry between 
the transverse and longitudinal directions for the optical 
potential in the lin_Llin lattice. In the physical picture 
of Sisyphus cooling p|]l4f| , cooling ends once an atom is 
trapped in a single potential well and hence the steady- 
state temperature is proportional to the potential depth. 
As already discussed in Sec. IV this is about twice as 



large in the transverse direction as in the longitudinal. 
Hence, the transverse temperature would be expected to 
be about twice as large as the longitudinal. However, cor- 
relations between these two directions tend to equilibrate 
the temperatures and in the simulations we therefore find 
T x to be only about 1.4 times larger than T z . 

We now turn to the jumping regime, |A| < 10r. We 
find a completely different behavior of the temperature. 
This is a consequence of the increasing contribution of 
absorption and spontaneous emission to the force fluctu- 
ations pSfl . Physically, an atom experiences many pho- 
ton recoils during one elementary cooling process and is 
thus more likely to escape from the trapping potential. 
Hence, the atom can reach a steady state temperature 
larger than the potential well depth. 



6030 
4000 
2000 





(a) - 








^^-Jf 300 - 



2.5 5 7.5 10 12.5 15 17.5 




In Fig. we plot the increase of the temperature com- 
pared to the oscillating regime as a function of the lattice 
periods A^. In excellent agreement with Eq. ( |34| ) we find 
that Ti is proportional to Aq and to Xf . 



VI. SPATIAL DIFFUSION OF THE ATOMIC 
CLOUD 

We now turn to the study of the spatial diffusion of the 
atomic cloud. In the simulations we calculate the average 
square position (rf (t)} in each direction over the whole 
cloud. 

In the thermalization phase the hot atoms follow al- 
most ballistic trajectories and the cloud expands rapidly. 
For longer times the expansion reaches a normal diffusion 
regime where 



(t)) = 2D Sl t + At 



(36) 



and Arf is a constant depending on the initial space and 
velocity distribution. Note that D$i does not depend on 
this initial distribution as wc verified in the simulations. 
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FIG. 8. Transverse (a) and longitudinal (b) spatial diffu- 
sion coefficients as a function of A for various A' (the values 
of [Ao[/av are plotted in the graphs) and for — 45°. The line 
corresponds to the analytical fit, Eq. (fj7|), for Aq = — 450o) r . 



FIG. 7. Variation of the (a) transverse and (b) longitudinal 
temperature as a function of the lattice spatial periods for 
small detunings. T° denotes the temperature in the domain 
of large A, A/T = —2 (triangles), A/T = —3 (squares), and 
|Anl/kV — 300, 600. The lines are linear fits of the form of 
Eq.H 



In Fig. g we plot the transverse and longitudinal spa- 
tial diffusion coefficients versus the lattice detuning for 
various light shifts and for a given lattice angle. Figure ^ 
clearly shows two domains where the spatial diffusion co- 
efficient behaves differently. For small detunings, D$i de- 
creases rapidly with |A| and increases with |A |, whereas 
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for large detunings, D$i increases with |A| and does not 
depend on A' except for A' = — 150w r . The latter cor- 
responds to a relatively shallow potential and the system 
is close to the transition to anomalous diffusion f^,[7|,^8| . 

Fitting the numerical results with an expression of the 
form of Eq. (|32|), we find: 



D$x = 0.50-^- f ^~ 

M \k x 

< — TXT + l-73^[1.16(Ufc) 2 + 0.83(fc z /fc) 2 ] 
ujr A 1 



On the contrary, in the oscillating regime we find: 



D 



Sz 



1.30 



M \k 



|A | r 



--- T 1.00V[«-83(fc z /fc) 2 + l.l6(k x /k) 2 ] 

wr a r 



(37) 



The coefficients of this fit are in good quantitative agree- 
ment with the analytical result ( |32| ) . The main difference 
is the factor of 1.39 which amounts to the difference be- 
tween the longitudinal and the transverse temperatures 
as found in Sec. [v|. Figure || shows that Eq. ( p7j) y ields 
an excellent fit to the numerical data. From Eq. (|1C|) and 
Eqs. (^2|), ( |37j ) we observe that the two domains where 
Dg behaves differently correspond to the jumping and 
the oscillating regimes. 

In the jumping regime we find: 



D Sx = 0.025 A 2 r , 
D Sz = 0.066 A 2 r , 



(38a) 
(38b) 



in good qu alit ative agreement with our physical discus- 
sion in Sec. 



Ill, see Eq. (13). In particular, D$i is propor- 
tional to A 2 as shown in Fig. ^ and proportional to the 
optical pumping rate T' . The different values of the nu- 
merical coefficients indicate that the trapping is stronger 
in the transverse direction than in the longitudinal, in 
agreement with our discussion in Sec. IV. 



D Sx = 



h 

M 



1+0.72 



Dsn = 1.08- 



1 + 1.39 I -- 

A x 



(39a) 



y. (39b) 



In this regime, D$ is proportional to A as expected 
from Eq. ( |l4| ) and the angular dependence is qualita- 
tively given by Eq. (^). In Fig. [l(j we show this an- 
gular dependence in the oscillating regime. Here, the 
transverse and longitudinal directions are not indepen- 
dent because the potential wells significantly deflect the 
trajectories of atoms travelling over many optical poten- 
tial wells. The dependence of D$i on 9 thus contains 
both \ x and X z . This is dramatically different to the 
situation in the jumping regime where atoms only jump 
between adjacent wells. Equation (|3^) also confirms that 
for A,; <C Aj, Dsi does not depend on Ai as predicted in 



Sec. III. 
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FIG. 10. Transverse (x) and longitudinal (z) spatial diffu- 
sion coefficients divided by A| as a function of 9 for A = — 30r 
and Aq = — 300cj r . The points correspond to the numerical 
simulations and the lines correspond to Eq. (^) (oscillating 
regime) . 
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FIG. 9. Transverse (x) and longitudinal (z) spatial diffu- 
sion coefficients divided by the corresponding square spatial 
period Xf as a function of 9 for A' — — 750w r and A = — 3r 
(jumping regime). 



The theoretical model of Sec. |iyj was based on the de- 
scription of the atomic dynamics by a Brownian motion 
model. Let us now further investigate the validity of 
such a description by testing in the numerical simula- 
tions some characteristics of Brownian motion. We par- 
ticularly perform a direct numerical measurement of the 
friction coefficients and test the validity of the Einstein 
relation (||). 

In order to probe the atomic dynamics we submit the 
atoms to a constant, space and velocity independent force 
in addition to the forces due to the atom-light interaction. 
In an experiment this could be provided simply by gravity 
or, for example, by the radiation pressure force of an 
additional weak laser beam. In a Brownian motion model 
such a constant force F^ c ^ = Fx + Fz C ^e z will give 
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rise to a constant mean velocity (v) 
atomic cloud with 



v x Vx + " 2 e z of the 



(40) 



Because of the linearity of the Brownian equations of 
motion the kinetic temperature and the spatial diffusion 
coefficients are not changed. 

Adding such a constant force in the numerical simula- 
tions along the i-direction, we observed in fact that the 
atomic cloud experiences a drift in this direction at a 
constant velocity. In the ideal case of a pure Brownian 
motion any amplitude of can be used, but in the 
case of Sisyphus cooling the friction coefficient is veloc- 

(c) 

ity dependent. In order to get Vi proportional to F> \ 
it is thus essential to use a small enough force which in- 
duces a global drift much smaller than the width of the 
velocity distribution. Under this condition the tempera- 
ture and the spatial diffusion do not depend on F^ c ' and 
Eq. (^) can be used to numerically find unique friction 
coefficients cti. 
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motion model. Note that the eigen-directions of the mo- 
tion are the x and z directions in good agreement with 
the theoretical model. 

In order to better understand this result let us briefly 
return to the model developed in Sec. |y| As has been 
shown in the previous section, the numerically obtained 
diffusion coefficient agrees well with the analytical result 
(|32|). In the derivation of the latter we have assumed 
that the cloud of atoms can be split into a trapped frac- 
tion and into a free fraction. Only the free atoms were 
taken into account for the spatial diffusion. Similar argu- 
ments must also be considered for the friction coefficient. 
Adding a small constant force will leave a trapped par- 
ticle in its initial potential well, and therefore the frac- 
tion of trapped atoms does not contribute to the mean 
velocity (v) of the cloud. Thus, the measured friction 
coefficient using Eq. ([lo]) should be given by the ana- 
lytical result (E4j) divided by the fraction of free atoms. 
Hence, for the measured values of the spatial diffusion 
and of the friction both sides of Eq. (||) are corrected 
by the same factor. In other words, the Einstein rela- 
tion holds because the measured quantities only involve 
the freely travelling atoms for which a Brownian motion 
model works well. 

Therefore, we obtain an analytic fit to the measured 
friction coefficient by inserting Eqs. (j3f]) and Q37|) into 
Eq. (§). We find 



a. 
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1 + 0.4 (fcf + 1.39*2) 
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(41b) 
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FIG. 11. Ratio of the friction coefficients calculated with 
Eq. dio| ) and with Eq. (|) for various A and for |Ao|/«r = 150, 
300, 450. 

We can then compare these numerical results with 
the friction coefficient obtained via the Einstein relation 
(|J) using the numerically found values of the tempera- 
ture (see Sec. |v|) and of the spatial diffusion coefficient 
(Sec. VI). In Fig. [ll] we plot the ratio of these two friction 
coefficients and find deviations of about 15%. Hence, the 
dynamics of an atomic cloud in an optical lattice is in 
reasonable agreement with a two-dimensional Brownian 



These approximate expressions are compared with the 
numerically obtained values of the friction coefficient in 
Fig. [l^. We see that there is excellent agreement both 
qualitatively and quantitatively. 

The behavior of on is again different in the jumping 
and the oscillating regime. In the jumping regime, a, is 
proportional to |A|/T and approximately independent of 



A §0: 



2hkt 



(42a) 
(42b) 



However, Fig. [12] exhibit a small dependence of on versus 
Aq in the jumping regime and this is not forecasted by 
the model. In fact, the kinetic temperature is not pro- 
portional but linear in Aq and this induce a dependence 
of the spatial distribution P in Aq. Coefficients Cj are 
thus Ag-dependent and this can explain the discrepancy. 

In the oscillating regime, oii is proportional to T' = 
TAq/A and depends on both X x and A 2 as discussed in 
Sec. Vi: 



10 



2M 



a z = 



l + 0.72fc2/fc2 A 

1.38M rA 
l + 1.39fc2/fc2~A" 



(43a) 
(43b) 



The expression found in the simulations is in good quali- 
tative agreement with the expression derived in the theo- 
retical model Eq. (fl3|). Note however that coefficients C; 
are not given by Eq. [30| but are to be calculated consid- 
ering the free atoms which contribute to spatial diffusion 
only. 




FIG. 12. Transverse (a) and longitudinal (b) friction coef- 
ficients versus A (data points) and the curves corresponding 
to Eq. (Eh for 9 = 45° and for A' /aj r = -150, -450, -750. 



VIII. CONCLUSIONS 

In conclusion, we performed a systematic study of the 
behavior of an atomic cloud in a lin_Llin optical lattice 
with the help of semi-classical Monte-Carlo simulations. 
We explored a broad range of lattice parameters includ- 
ing the jumping and the oscillating regime. 

The temperature was found to be linear in the poten- 
tial depth and independent of the laser detuning and of 
the lattice angle in a broad range of parameters. We have 
shown that the temperature is anisotropic, the transverse 
one being larger than the longitudinal one by a factor of 
1.4. All these results are well explained with the help of 
the physical picture of Sisyphus cooling and are in good 
agreement with experimental results. 

The spatial difusion D$ was studied in the regime of 
normal diffusion. The behavior of D$ differs significally 



in the jumping and in the oscillation regimes. In the first, 
Ds decreases with |A| and increases with |A |. In the 
second, Ds increases linearly with |A| and does not de- 
pend on |Ag|. The behavior of Ds as a function of the 
lattice spatial periods is also different in both regimes: 
whereas Dsi is proportional to Af in the jumping regime, 
it is a function of A^ as well as X z in the oscillating regime. 
This reveals correlations between the transverse and lon- 
gitudinal directions of the lattice. 

By adding a constant force in the Monte-Carlo sim- 
ulations we could numerically measure the friction co- 
efficients and we showed that the Einstein relations are 
fulfilled. This supports a description of the dynamics in 
terms of Brownian motion. The friction coefficient a, is 
proportional to A and 1/Af in the jumping regime. In 
the oscillating regime, cti is proportional to A' and 1/A, 
and the dependence on the lattice geometry involves both 
l/X x and 1/A Z . 

The numerical results have been found to be in good 
agreement with a simple theoretical model based on a 
semi-classical approach. We derive the steady-state tem- 
perature, the friction force, and the spatial diffusion from 
a model of Brownian motion taking into account atomic 
localization in the optical potential wells. To explain 
the measured friction and spatial diffusion, the atomic 
cloud must be split into a trapped part and a free part. 
While in general both parts contribute to the internal 
and external dynamics, only the free fraction of atoms is 
responsible for the observed expansion of the cloud and 
for the drift of the center of mass under the influence of 
a weak constant force. 

The spatial diffusion of atomic clouds in optical lat- 
tices is usually studied in pump-probe spectroscopy ex- 
periments, using the properties of the Rayleigh line [j30| . 
However, the validity of this method has never been 
proven. We expect that the models discussed here will 
serve to this verification by providing a systematic the- 
oretical study of the directly measured spatial diffusion 
coefficients. 

Finally, it should be noted that our restriction to 
Ox = y has been guided by experimental restrictions 
but is not necessary. Indeed, when this condition is not 
fulfilled, the unbalanced radiation pressure induces a fast 
escape of the atomic cloud from the optical lattice which 
makes experimental investigations difficult. Neverthe- 
less, this situation could be of great interest because it 
offers the opportunity of studying optical lattices with 
different spatial periods along the x and y axes and in- 
teresting anisotropic effects could be found. 
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